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Abstract 

While the seismic quality factor and phase lag are defined solely by the bulk proper- 
ties of the mantle, their tidal counterparts are determined both by the bulk properties 
and the size effect (self-gravitation of a body as a whole). For a qualitative estimate, we 
model the body with a homogeneous sphere, and express the tidal phase lag through 
the lag in a sample of material. Although simplistic, our model is sufficient to under- 
stand that the lags are not identical. The difference emerges because self-gravitation 
pulls the tidal bulge down. At low frequencies, this reduces strain and the damping 
rate, and makes tidal damping less efficient in larger objects. At higher frequencies, 
competition between self-gravitation and rheology becomes more complex, though for 
sufficiently large superearths the same rule applies: the larger the planet, the weaker 
tidal dissipation in it. Being negligible for small terrestrial planets and moons, the 
difference between the seismic and tidal lagging (and likewise between the seismic and 
tidal damping) becomes very considerable for large exoplanets (superearths). In those, 
it is much lower than what one might expect from using a seismic quality factor. 

The tidal damping rate deviates from the seismic damping rate especially in the 
zero- frequency limit, and this difference takes place for bodies of any size. So the 
equal in magnitude but opposite in sign tidal torques, exerted on one another by the 
primary and the secondary, have their orbital averages going smoothly through zero as 
the secondary crosses the synchronous orbit. 

We describe the mantle rheology with the Andrade model, allowing it to lean toward 
the Maxwell model at the lowest frequencies. To implement this additional flexibility, 
we reformulate the Andrade model by endowing it with a free parameter £ which is the 
ratio of the Andrade timescale to the viscoelastic Maxwell time of the mantle. Some 
uncertainty in this parameter's frequency dependence does not influence our principal 
conclusions. 



1 The goal and the plan 

As the research on exoplanetary systems is gaining momentum, more and more accurate 
theoretical tools of planetary dynamics come into demand. Among those tools are the 
methods of calculation of tidal evolution of both orbital and rotational motion of planets 



1 



and their moons. Such calculations involve two kind of integral parameters of celestial bodies 
- the Love numbers and the tidal quality factors. The values of these parameters depend 
upon the rheology of a body, as well as its size, temperature, and the tidal frequency. 

It has recently become almost conventional in the literature to assume that the tidal 
quality factor of superearths should be of the order of one hundred to several hundred 
(Carter et al. 2011, Leger et al. 2009). Although an acceptable estimate for the seismic Q , 
this range of numbers turns out to fall short, sometimes by orders of magnitude, of the tidal 
Q of superearths. 

In our paper, the frequency dependence of tidal damping in a near-spherical homogeneous 
body is juxtaposed with the frequency dependence of damping in a sample of the material 
of which the body consists. For brevity, damping in a sample will be termed (somewhat 
broadly) as "seismic damping". 

We shall demonstrate that, while the tidal Q of the solid Earth happens not to deviate 
much from the solid-Earth seismic Q , the situation with larger telluric bodies is considerably 
different. The difference stems from the presence of self-gravitation, which suppresses the 
tidal bulge and thereby acts as extra rigidity - a long-known circumstance often neglected in 
astronomical studiesjj Due to self-gravitation ("size effect"), tidal damping in superearths 
is much less efficient than in earths, and the difference may come to orders of magnitude, as 
will be demonstrated below. Thus, while the seismic Q of a superearth may be comparable 
to the seismic Q of the solid Earth, the tidal Q of a superearth may exceed this superearth's 
seismic Q greatly. This is the reason why it is inappropriate to approximate superearths' 
tidal quality factors with that of the solid Earth. 

We shall show that the difference between the frequency dependence of the tidal Q factor 
and that of the seismic Q may explain the "improper" frequency dependence of the tidal 
dissipation rate measured by the lunar laser ranging (LLR) method. We also shall point 
out that the correct frequency dependence of the tidal dissipation rate, especially at low 
frequencies, plays an important role in modeling the process of entrapment into spin-orbit 
resonances. In greater detail, the latter circumstance will be discussed in Efroimsky (2012). 

The rate of the "seismic damping" (a term that we employ to denote also damping in a 
sample of the material) is defined, at each frequency, by the material's rheology only, i.e., by 
the constitutive equation linking the strain and stress at this frequency. The rate of the tidal 
damping however is determined both by the rheology and by the intensity of self-gravitation 
of the body. At a qualitative level, this can be illustrated by the presence of two terms, 1 
and 19/x(oo)/(2pgi?) , in the denominator of the expression for the static Love number &2 
of a homogeneous sphere. Here /U.(oo) denotes the relaxed shear modulus, g signifies the 
surface gravity, while p and R stand for the mean density and the radius of the body. The 
first of these terms, 1 , is responsible for the size effect (self-gravitation), the second for the 
bulk properties of the medium. Within the applicability realm of an important theorem 
called elastic-viscoelastic analogy (also referred to as the correspondence principle), the same 
expression interconnects the Fourier component k 2 (x) °f the time derivative of the Love 
number with the Fourier component p>(x) °f the stress- relaxation function at frequency x ■ 

1 Including the size effect via ki is common. Unfortunately, it is commonly assumed sufficient. This 
treatment however is inconsistent in that it ignores the inseparable connection between the Love number 
and the tidal quality factor (or the tidal phase lag). In reality, both the Love number and the sine of the 
tidal lag should be derived from the rheology and geometry of the celestial body, and cannot be adjusted 
separately from one another. 



2 



This renders the frequency-dependence of the tangent of the tidal lag, which is the negative 
ratio of the imaginary and real parts of ^(x) • 

This preliminary consideration illustrates the way rheology enters the picture. First, the 
constitutive equation defines the frequency dependence of the complex compliance, J(x) , 
and of the complex rigidity JJl{x) = 1/J(x) ■ The functional form of this dependence deter- 
mines the frequency dependence of the complex Love number, ki(x) ■ The latter furnishes 
the frequency dependence of the products |&z(x)| sinez(x) which enter the tidal theory. 

In Section we briefly recall the standard description of stress-strain relaxation and 
dissipation in linear media. In Section [3], we describe a rheological model, which has proven 
to be adequate to the experimental data on the mantle minerals and partial melts. The 
goal of the subsequent sections will be to build the rheology into the theory of bodily tides, 
and to compare a tidal response of a near-spherical body to a seismic response rendered by 
the medium. Finally, several examples will be provided. Among these, will be the case of 
the Moon, whose "improper" tidal-dissipation frequency-dependence finds an explanation 
as soon as the difference between the seismic and tidal friction is brought to light. In the 
closing section, we shall compare our results with those obtained by Goldreich (1963). 

2 Formalism 

Everywhere in this paper we shall take into consideration only the deviatoric stresses and 
strains, thus neglecting compressibility. 

2.1 Compliance and rigidity. 

The standard linear formalism in the time domain 

The value of strain in a material depends only on the present and past values taken by the 
stress and not on the current rate of change of the stress. Hence the compliance operator J 
mapping the stress cr 7 „ to the strain w 7 „ must be just an integral operator, linear at small 
deformations: 

2u lv (t) = J(t) <v = / J(t - t') a lu {t')dt' , (1) 

J — oo 

where t' <t , while overdot denotes d/dt' . The kernel J(t — t') is termed the compliance 
function or the creep-response function. 
Integration by parts renders: 

2u lu {t) = J(t)a lv = J(0) <v(t) - J(oo) <v(-oo) + f j (t - t') <v(t') dt' . (2) 

J — oo 

As the load in the infinite past may be set zero, the term containing the relaxed compliance 
J(oo) may be dropped. The unrelaxed compliance J(0) can be absorbed into the integral 
if we agree that the elastic contribution enters the compliance function not as 

J(t — t') = J(0) + viscous and hereditary terms , (3) 

but as 

J(t — t') = J(0) 0(t — t') + viscous and hereditary terms . (4) 
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The Heaviside step-function Q(t — t ') is set unity for t — t' > , and zero for t — t' < , 
so its derivative is the delta-function S(t — t') . Keeping this in mind, we reshape fl2]) into 



2ury U (t) = J(t) o> = / J (t-t') a^(t')dt' , with J(t-t') containing J(0)6(t-t') . (5) 

Inverse to the compliance operator 

2u lu = J a lv . (6) 

is the rigidity operator 

a yu = 2pi u lu . (7) 

In the presence of viscosity, operator jl is not integral but is integrodifferential, and thus 
cannot be expressed as a JU (t) = 2 fi(t — t') u JU (t') dt' . It can though be written as 

t 



a lu {t) = 2 / fi(t - t')u lv {t')dt' , (8) 

J — oo 

if its kernel, the stress-relaxation function fi(t — t ') , is imparted with a term 2 1] 5(t — 
t ') , integration whereof renders the viscous portion of stress, 2 77 u lv . The kernel also 
incorporates an unrelaxed part /x(0) Q(t — t ') , whose integration furnishes the elastic portion 
of the stress. The unrelaxed rigidity /i(0) is inverse to the unrelaxed compliance J(0) . 

Each term in fi(t — t') , which is neither constant nor proportional to a delta function, 
is responsible for hereditary reaction. 

2.2 In the frequency domain 

To Fourier-expand a real function, nonnegative frequencies are sufficient. Thus we write: 

f°° ixt f°° ixt 

= / cv(x) e dx and u lv (t) = / w 7 „(x) e d\ , (9) 
Jo Jo 

where the complex amplitudes are 

MX) = ^(X) e^ CT(x) , MX) = e^ (x) , (10) 

while the initial phases <f a (x) an d <Pu(x) are se t to render the real amplitudes ov(x„) 
and M 7 j/(x n ) non- negative. To ensure convergence, the frequency is, whenever necessary, 
assumed to approach the real axis from below: Im(x) — > 0— . 

With the same caveats, the complex compliance J(x) is introduced as the Fourier image 
of the time derivative of the creep-response function: 



J(x)e lXT d X =J(r) . (11) 

Jo 

The inverse expression, 

/>oo 

AX) = / j(r)e- lXT dr , (12) 
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is often written down as 



J( X ) = J(0) + i X [ J(r) - J(0)6(r) ] dr . (13) 

Jo 

For causality reasons, the integration over r spans the interval [ 0, oo ) only. Alternatively, 
we can accept the convention that each term in the creep-response function is accompanied 
with the Heaviside step function. 

Insertion of the Fourier integrals f)9l- [TTj) into (JT|) leads us to 

oo • , p oo • , 



2 / ?v(x) e dx = <v(x) Ax) e rfx , (14) 
Jo Jo 

whence we obtain: 

2 = -Ax) ^(x) • (15) 

Expressing the complex compliance as 

J(X) = \J(X)\ exp [-»%)] , (16) 

where 

, Im | J (x) 1 , . 

tan 5 x = 1- J , 17 

K J Ke[J{ X )] 1 ^ 

we see that 5(x) is the phase lag of a strain harmonic mode relative to the appropriate 
harmonic mode of the stress: 

<pM = Mx) - S( X ) . (is) 

2.3 The quality factor(s) 

In the linear approximation, at each frequency x t ne average (per period) energy dissipation 

rate (E (x)) is defined by the deformation at that frequency only, and bears no dependence 
upon the other frequencies: 

(E(X)} = -^^ (19) 

or, the same: 

AE cyde ( X ) = - , (20) 

AE cyc i e (x) being the one-cycle energy loss, and Q(x) being the quality factor related to the 
phase lag at the frequency x ■ It should be clarified right away, to which of the lags we are 
linking the quality factor. When we are talking about a sample of material, this lag is simply 
5(x) introduced above as the negative argument of the appropriate Fourier component of the 
complex compliance - see formulae (TTTl- fTBl . However, whenever we address tide, the quality 
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factor becomes linked (via the same formulae) to the tidal phase lag e(x) . Within the same 
rheological model, the expression for e(x) differs from that for 5(x) , because the tidal lag 
depends not only upon the local properties of the material, but also upon self-gravitation of 
the body as a whole. 

The aforementioned "seismic-or-tidal" ambiguity in definition of Q becomes curable as 
soon as one points out to which kind of deformation the quality factor pertains. More serious 
is the ambiguity stemming from the freedom in defining E pea k(x) ■ 

If E pea k{x) i n CGS-EU]) signifies the peak energy stored at frequency %■> the resulting 
quality factor is related to the lag via 

Q- 1 = sin 151 (21) 

^ energy 11 y ' 

(not tan \8\ as commonly believed - see the calculation in the Appendix to Efroimsky 2012). 

If however E pea k(x) is introduced as the absolute maximum of work carried out on the 
sample at frequency x over a time interval through which the power stays positive, then 
the appropriate Q factor is connected to the lag via 

^ i tan 151 , , 

Q~ \ = t > 22 

^ work „ / — . . \ . . . ' 



\5\ tan |^| 



7T 
, 1 

as was shown in 

The third definition of the quality factor (offered by Goldreich 1963) is 

QgoI^c^^W ■ ( 23 ) 

This definition, though, corresponds neither to the peak work nor to the peak energy. 
In the limit of weak lagging, all three definitions entail 

Q- 1 = \S\ + 0(5 2 ) . (24) 

For the lag approaching tt/2, the quality factor defined as ( J2TT) assumes its minimal value, 
Q energy = ^ ' wnne definition ( l22l) renders Q work = . The latter is natural, since in the 
considered limit the work performed on the system is negative, its absolute maximum being 
zerojf] 

In seismic studies or in exploration of attenuation in small samples, one's choice among 
the three definitions of Q is a matter of personal taste, for the quality factor is large and 
the definitions virtually coincide. 

In the theory of tides, the situation is different, because at times one has to deal with 
situations where the definitions of Q disagree noticeably - this happens when dissipation 
is intensive and Q is of order unity. To make a choice, recall that the actual quantities 



2 In Efroimsky & Williams (2009), E pea k(x) was miscalled "peak energy". However the calculation of Q 
was performed there for E pea k (x) introduced as the peak work. 

3 As Q < 2ir implies E pea k < AE , such small values of Q are unattainable in the case of damped free 
oscillations. Still, Q can assume such values under excitation, tides being the case. 
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entering the Fourier expansion of tides over the modes to lmpq are the products^ 

h sine, = ^(w, mM ) Bine { (w frnM ) , (25) 

where ^(w Zmp(? ) are the dynamical analogues to the Love numbers. It is these products that 
show up in the Impq terms of the expansion for the tidal potential (force, torque). From 
this point of view, a definition like ( )2T1) would be preferable, though this time with the tidal 
lag e instead of the seismic lag 5 : 



Q 1 = sin | e t | (26a) 



or, in a more detailed manner: 



Qi'^imJ = sin 1 6,(0^)1 . (26b) 

Under this definition, one is free to substitute k\ sine z with ki/Q^ The subscript / 
accompanying the tidal quality factor will then serve as a reminder of the distinction between 
the tidal quality factor and its seismic counterpart. 

While the notion of the tidal quality factor has some illustrative power and may be 
employed for rough estimates, calculations involving bodily tides should be based not on the 
knowledge of the quality factor but on the knowledge of the overall frequency dependence 
of products ki sine ; = &/(c^ mp(J ) sine^w^^) . Relying on these functions would spare one 
of the ambiguity in definition of Q and would also enable one to take into account the 
frequency dependence of the dynamical Love numbers. 



3 The Andrade model and its reparameterisation 

In the low-frequency limit, the mantle's behaviour is unlikely to differ much from that of the 
Maxwell body, because over timescales much longer than 1 yr viscosity dominates (Karato & 
Spetzler 1990). At the same time, the accumulated geophysical, seismological, and geodetic 
observations suggest that at shorter timescales inelasticity takes over and the mantle is 
described by the Andrade model. However, the near-Maxwell behavior expected at low 
frequencies can be fit into the Andrade formalism, as we shall explain below. 



3.1 Experimental data: the power scaling law 

Dissipation in solids may be effectively modeled using the empirical scaling law 

sind = (S X f P , (27) 

£ being a constant having the dimensions of time. This "constant" may itself bear a (typ- 
ically, much slower) dependence upon the frequency x • The dependence of £ on the tem- 
perature is given by the Arrhenius law (Karato 2008). 

4 A historical tradition (originating from Kaula 1964) prescribes to denote the tidal phase lags with 
e impq ' wnne keeping for the dynamical Love numbers the same notation as for their static predecessors: ki . 
These conventions are in conflict because the product fc; sin e tolp(? is the negative imaginary part of the 
complex Love number ki . More logical is to use the unified notation as in (f25j) . At the same time, it should 
not be forgotten that for triaxial bodies the functional form of the dependence of fc; on frequency is defined 
not only by I but also by m, p, q . In those situations, one has to deal with k t s ^ ne imvo > see Section |U 
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Experiments demonstrate that the power dependence (1271) is surprisingly universal, with 
the exponential p robustly taking values within the interval from 0.14 to 0.4 (more often, 
from 0.14 to 0.3). 

For the first time, dependence ( 12 7p was measured on metals in a lab. This was done by 
Andrade (1910), who also tried to pick up an expression for the compliance compatible with 
this scaling law. Later studies have demonstrated that this law works equally well, and with 
similar values of p, both for silicate rocks (Weertman & Weertman 1975, Tan et al. 1997) 
and ices (Castillo- Rogez 2009, McCarthy et al 2007). 

Independently from the studies of samples in the lab, the scaling behaviour ( 1271) was 
obtained via measurements of dissipation of seismic waves in the Earth (Mitchell 1995, 
Stachnik et al. 2004, Shito et al. 2004). 

The third source of confirmation of the power scaling law came from geodetic experiments 
that included: (a) satellite laser ranging (SLR) of tidal variations in the J2 component of 
the gravity field of the Earth; (b) space-based observations of tidal variations in the Earth's 
rotation rate; and (c) space-based measurements of the Chandler wobble period and damping 
(Benjamin et al. 2006, Eanes & Bettadpur 1996, Eanes 1995)0 

While samples of most minerals furnish the values of p lying within the interval 0.15 — 
0.4, the geodetic measurements give 0.14 — 0.2. At least a fraction of this difference may 
be attributed to the presence of partial melt, which is known to have lower values of p 
(Fontaine et al. 2005). 

On all these grounds, it is believed that mantles of terrestrial planets are adequately 
described by the Andrade model, at least in the higher frequency band where inelastic- 
ity dominates (Gribb & Cooper 1998, Birger 2007, Efroimsky & Lainey 2007, Zharkov & 
Gudkova 2009). Some of the other models were considered by Henning et al. (2009). 

The Andrade model is equally well applicable to celestial bodies with ice mantles (for 
application to Iapetus see Castillo- Rogez et al. 2011) and to bodies with considerable hy- 
dration in a silicate mantled The model can also be employed for modeling of the tidal 
response of the solid parts of objects with significant liquid-water layersE) 

3.2 The Andrade model in the time domain 

The compliance function of the Andrade body (Cottrell & Aytekin 1947, Duval 1978), 

J(t-t') = [ J + (t-t') a (5+ it-t') r]- 1 ] e(t-t') , (28) 

contains empirical parameters a and /3 , the steady-state viscosity rj , and the unrelaxed 
compliance J = J(0) = l//i(0) = 1/fi. We endow the right-hand side of (|28|) with the 
Heaviside step-function O (£ — £') , to ensure that insertion of ( )28|) into ([5]), with the sub- 
sequent differentiation, yields the elastic term J 5{t — t') under the integral. The model 

5 It should be noted that in reality the geodetic measurements were confirming the power law ([27)) not 
for the seismic lag 8 but for the tidal lag e , an important detail to be addressed shortly. 

6 Damping mechanisms in a wet planet will be the same as in a dry one, except that their efficiency will 
be increased. So the dissipation rate will have a similar frequency dependence but higher magnitude. 

7 In the absence of internal oceans, a rough estimate of the tidal response can be obtained through model- 
ing the body with a homogeneous sphere. However the presence of such oceans makes it absolutely necessary 
to calculate the overall response through integration over the solid and liquid layers. As demonstrated by 
Tyler (2009), tidal dissipation in internal ocean layers can play a big role in rotational dynamics of the body. 
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allows for description of dissipation mechanisms over a continuum of frequencies, which is 
useful for complex materials with a range of grain sizes. 

The Andrade model can be thought of as the Maxwell model equipped with an extra 
term (t — t') a (3 describing hereditary reaction of strain to stress. The Maxwell model 



(Maxwell) 
J \t~t') 



[j + {t-t') rr 1 ] e(*-o 



(29) 



is a simple rheology, which has a long history of application to planetary problems, but 
generally has too strong a frequency dependence at higher frequencies where inelasticity 
becomes more efficient than viscosity (Karato 2008). Insertion of (T2"9"j) into (j3J) renders strain 
consisting of two separate inputs. The one proportional to J implements the instantaneous 
(elastic) reaction, while the one containing 7/ -1 is responsible for the viscous part of the 
reaction. 

Just as the viscous term (t — t') rj^ 1 showing up in f[2"8l - 12"§|) . so the inelastic term 
(t — t f ) a (3 emerging in the Andrade model (|28|) is delayed - both terms reflect how the 
past stressing is influencing the present deformation. At the same time, inelastic reaction 
differs from viscosity both mathematically and physically, because it is produced by different 
physical mechanisms. 

A disadvantage of the formulation (1281) of the Andrade model is that it contains a pa- 
rameter of fractional dimensions, /3 . To avoid fractional dimensions, we shall express this 
parameter, following Efroimsky (2012), as 



P 



A 



(30a) 



the new parameter r A having dimensions of time. This is the timescale associated with the 
Andrade creep, wherefore it may be named as the 'Andrade time" or the "inelastic time". 
Another option is to express (3 as 



a , , — 1 a 

M 



where the dimensionless parameter ( is related through 



7\ 



c = — 



(30b) 



(31) 



to the Andrade timescale r A and to the Maxwell time 



t, 



V T 

— = rj J . 



(32) 



In terms of the so-introduced parameters, the compliance assumes the form of 



J(i-t') 



J 



t-t' 



l \ a 



+ 



t-t' 



G(t-t') 



(33a) 



J 



t-t' 

C r M 



t-t' 



Q(t-t') 



(33b) 
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For r 4 <C t m (or, equivalently, for ( <C 1 ), inelasticity plays a more important role than 
viscosity. On the other hand, a large r A (or large ( ) would imply suppression of inelastic 
creep, compared to viscosity. 

It has been demonstrated by Castillo- Rogez that under low stressing (i.e., when the 
grain-boundary diffusion is the dominant damping mechanism - like in Iapetus) obeys 
the relation 

(3 « J T ~ a = J l ~ a rT a = r]~ a , (34a) 

(see, e.g., Castillo- Rogez et al. 2011). Comparing this to fl30|) . we can say that the Andrade 
and Maxwell timescales are close to one another: 

r A « r M (34b) 

or, equivalently, that the dimensionless parameter ( is close to unity: 

C « 1 . (34c) 

Generally, we have no reason to expect the Andrade and Maxwell timescales to coincide, 
nor even to be comparable under all possible circumstances. While (1341) may work when 
inelastic friction is determined mainly by the grain-boundary diffusion, we are also aware of 
a case when the timescales r A and r M differ considerably. This is a situation when stressing 
is stronger, and the inelastic part of dissipation is defined mainly by dislocations unpinning 
(i.e., by the Andrade creep). This is what happens in mantles of earths and superearths. 

On theoretical grounds, Karato & Spetzler (1990) point out that the dislocation- unpinning 
mechanism remains effective in the Earth's mantle down to the frequency threshold Xo ~ 
1 yr -1 . At lower frequencies, this mechanism becomes less efficient, giving way to viscosity. 
Thus at low frequencies the mantle's behaviour becomes closer to that of the Maxwell bodyj!l 
This important example tells us that the Andrade time r A and the dimensionless parame- 
ter ( may, at times, be more sensitive to the frequency than the Maxwell time would be. 
Whether r A and ( demonstrate this sensitivity or not - may, in its turn, depend upon the 
intensity of loading, i.e., upon the damping mechanisms involved. 

3.3 The Andrade model in the frequency domain 

Through (TT2]) . it can be demonstrated (Findley et al. 1976) that in the frequency domain 
the compliance of an Andrade material reads as 

J(X) = J + P (ix)~ a r (1 + a) — (35a) 

vx 

= J [1 + (i X T A r a T(l + a) - i (xrj- 1 } , (35b) 
= J [1 + (i X (T M r a r(l + a) - i{ X r M y l ] , (35c) 



8 Using the Andrade model as a fit to the experimentally observed scaling law (j!?T)) . we see that the 
exponential p coincides with the Andrade parameter a < 1 at frequencies above the said threshold, and 
that p becomes closer to unity below the threshold - see subsection 13.41 below. 
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X being the frequency, and T denoting the Gamma function. The imaginary and real parts 
of the complex compliance are: 

Tm[J(x)} = ~ — ~ X~ a P sin ( ^) T(a + 1) (36a) 



rjx ^ 2 

- JixrJ' 1 - J(xr A r a sm(^) T(a + 1) (36b) 

-J{XT M r l - J(x(r M y a sm(^) T(a + 1) (36c) 



and 

ne[J{x)] = J + x~ Q /3 cos(^) I> + 1) (37a) 
= J + J( X r A )- a c os(^) r(a + 1) (37b) 

= J + J(xCr M r a cos(^) r(a + 1) . (37c) 
The ensuing frequency dependence of the phase lag will look as0 

Xm[J{x)} ^x)' 1 + x""/? sin(«^)r(« + 1) 



^ e [J(X)] n _1 + x~° P cos fe) r (a + 1) 

f (yr:-)" r (a + 1) sin (^) z" 1 C + sin (^) r (a + 1) 



1+ f-^L-)" r(a + 1) cos(^) 1 + z~ a cos(^) r(a + 1) 

with z being the dimensionless frequency defined through 

z = X r A = X r M C • (39) 
Evidently, for (3 — > (that is, for ( — > oo or r A — > oo), expression (155]) approaches 

tan<5(x) = (r A/ x) _1 , (40) 
which is the frequency dependence appropriate to the Maxwell body. 

3.4 Low frequencies: from Andrade toward Maxwell 

The Andrade body demonstrates the so-called "elbow dependence" of the dissipation rate 
upon frequency. At high frequencies, the lag 5 satisfies the power law 

tan 5 ~ X ~ v , (41) 



9 In some publications (e.g., Nimmo 2008), formula (|38a|l is given as an expression for the inverse quality 
factor. This is legitimate when the latter is defined through ([23]) . 
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the exponential being expressed via an empirical parameter a , where < a < 1 for most 
materials. It follows from ( )38|) that at higher frequencies p = a , while at low frequencies 
p = 1 — a . 

The statement by Karato & Spetzler (1990), that the mantle's behaviour at low frequen- 
cies should lean toward that of the Maxwell body, can be fit into the Andrade formalism, 
if we agree that at low frequencies either a approaches zero (so p approaches unity) or £ 
becomes large (so r A becomes much larger than r M ). The latter option is more physi- 
cal, because the increase of r A would reflect the slowing-down of the unpinning mechanism 
studied in Ibid. 

One way or another, the so-parameterised Andrade model is fit to embrace the result from 
Ibid. Hence our treatment will permit us to describe both the high-frequency range where the 
traditional Andrade model is applicable, and the low-frequency band where the behaviour 
of the mantle deviates from the Andrade model toward the Maxwell body. Comparison of 
the two models in the frequency domain is presented in Figure 1. 
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Figure 1: Andrade and Maxwell models in the frequency domain. The plot shows the decadic 
logarithm of tan 5 , as a function of the decadic logarithm of the forcing frequency x (i n cycles 
s ). For the Andrade body, the tangent of phase lag is given by (l38b|) . with a = 0.2 and 
t a = t m = 10 10 s. For the Maxwell body, the tangent is rendered by (T40|) . with t m = 10 10 s. 
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4 Expanding a tidal potential or torque — over the tidal 
modes or over the forcing frequencies? 

Consider a binary system with mean motion n , and suppose that tidal dissipation in one of 
the bodies much exceeds that in its companion. Then the former body may be treated as 
the tidally perturbed primary, the latter being its tide-raising secondary. The sidereal angle 

and the spin rate of the primary will be expressed with 9 and 9, while the node, pericentre, 
and mean anomaly of the secondary, as seen from the primary] 10 ! will be denoted by f2 , co , 
and M. 

In the Darwin-Kaula theory, bodily tides are expanded over the modes 

u lmpq = (I - 2p) Cj + (I - 2p + q) M + m (f2 - 9) « (I - 2p + q) n - m 9 , (42) 

with /, m, p, q being integers. Dependent upon the values of the mean motion, spin rate, 
and the indices, the tidal modes 0Ji mpq may be positive or negative or zero. 

In the expansion of the tidal potential or torque or force, summation over the integer 
indices goes as Ylu=2 X^m=o X^=o - oo ■ For example, the secular polar component of 
the tidal torque will read as: 

oo I oo oo 

T = Yl k l( U lm Pq ) ^ e Mm Pq ) > ( 4 3) 

1=2 m=0 p=0 q= — oo 

where the ellipsis denotes a function of the primary's radius and the secondary's orbital 
elements. The functions ) are the dynamical analogues of the static Love numbers, 

while the phase lags corresponding to the tidal modes u lmpq are given by 

€ Mmpg) = ^Impg At l(hmpg\) = I ^Impg I At l ( Hmpg I ) S S n ^Impg = Xlmpg At l(Xlm pq ) Sgn W 

Here the positively defined quantities 

Xlmpg — I U lmpg I (45) 

are the forcing frequencies in the material, while the positively defined time lags At ; are 
their functions. 

Following Kaula (1964), the phase and time lags are often denoted with ei mpq and Ati mpq . 
For near-spherical bodies, though, the notations ei{xi mp g) an d Ati(xim Pq ) would be prefer- 
able, because for such bodies the functional form of the dependency e; mpg (x) is defined by 
Is, but is ignorant of the values of the other three indicesO The same therefore applies to 
the time lag. 

10 When the role of the primary is played by a planet and the role of the perturbing secondary is played 
by the host star, the argument of the pericentre of the star as seen from the planet, u> , differs by 7r from 
the argument of the pericentre of the planet as seen from the star. 

Also mind that in equation (|42j) the letter uj with the subscript Impq denotes, as ever, a tidal mode, 
while the same letter without a subscript stands for the periapse. The latter use of this letter in equation 
(|42[) is exceptional, in that elsewhere in the paper the letter ui , with or without a subscript, always denotes 
a tidal mode. 

11 Within the applicability realm of the elastic- viscoelastic analogy employed in subsection 15 . 21 below, the 
functional form of the complex Love number ki (\) of a near-spherical object is determined by index / solely, 
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The forcing frequencies in the material of the primary, Xi mp q > are positively defined. 
While the general formula for a Fourier expansion of a field includes integration (or summa- 
tion) over both positive and negative frequencies, it is easy to demonstrate that in the case 
of real fields it is sufficient to expand over positive frequencies only. The condition of the 
field being real requires that the real part of a Fourier term at a negative frequency is equal 
to the real part of the term at an opposite, positive frequency. Hence one can get rid of the 
terms with negative frequencies, at the cost of doubling the appropriate terms with positive 
frequencies. (The convention is that the field is the real part of a complex expression.) 

The tidal theory is a rare exception to this rule: here, a contribution of a Fourier mode 
into the potential is not completely equivalent to the contribution of the mode of an opposite 
sign. The reason for this is that the tidal theory is developed to render expressions for 
tidal forces and torques, and the sign of the tidal mode uJ t shows up explicitly in those 
expression. This happens because the phase lag in (1431) is the product ( 144)) of the tidal mode 
u im Pq and tlie positively defined time lag At lmpg . 

This way, if we choose to expand tide over the positively defined frequencies x ° n ly, we 
shall have to insert "by hand" the multipliers 



S ^ nUJ lm Pq = S § n 



(I - 2p + q) n - m9 (46) 



into the expressions for the tidal torque and force, a result to be employed below in formula 
f )65|) . The topic is explained in greater detail in Efroimsky (2012). 



5 Complex Love numbers and the elastic-viscoelastic 
analogy 

Let us recall briefly the switch from the stationary Love numbers to their dynamical counter- 
parts, the Love operators. The method was pioneered, probably, by Zahn (1966) who applied 
it to a purely viscous planet. The method works likewise for an arbitrary linear rheologi- 
cal model, insofar as the elastic-viscoelastic analogy (also referred to as the correspondence 
principle) remains in force. 

5.1 From the Love numbers to the Love operators 

A homogeneous near-spherical incompressible primary alters its shape and potential, when 

— * 

influenced by a static point-like secondary. At a point R = (R, A, 0), the potential due to a 
tide-generating secondary of mass M* ec , located at r* = (r*, A*, 0*) with r* > R, can be 

while the integers m, p, q show up through the value of the frequency: fc;(x) = ki(ximpq) ■ This applies 
to the lag too, since the latter is related to ki via (|62|) . 

For triaxial bodies, the functional forms of the frequency dependencies of the Love numbers and phase lags 
do depend upon m, p, q , because of coupling between spherical harmonics. In those situations, notations 
kimpq and ei m pq become necessary (Dchant 1987a, b; Smith 1974). The Love numbers of a slightly non- 
spherical primary differ from the Love numbers of the spherical reference body by a term of the order of the 
flattening, so a small non-sphericity can usually be ignored. 
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expressed through the Legendre polynomials P;(cos7) or the Legendre functions P; m (sin0): 

i 

00 C M* 00 / R \ 

w(r, ?*) = wiifi, n = - E 4 ^( cos 7) 

1=2 1=2 ^ ' 

00 ' z 

_GM^y^(R\ V{y|(2-yP to ( S ^)P ln ( S mf) cosm(A-A*) , (47) 

1=2 v 7 m=0 v ' 

G being Newton's gravitational constant, and 7 being the angle between the vectors r* 

— * 

and -R originating from the primary's centre. The latitudes 4>, <p* are reckoned from the 
primary's equator, while the longitudes A, A* are reckoned from a fixed meridian on the 
primary. 

The I th spherical harmonic Ui(r) of the resulting change of the primary's potential at 
an exterior point r is connected to the I th spherical harmonic Wi(R, r*) of the perturbing 
exterior potential via Ui{r) = (R/r) l+1 k[Wi(R, r*) , so the total change in the exterior 
potential of the primary becomes: 

u{r) = J2 U ^ = E ( - ) k i w i&> o • ( 48 ) 

1=2 1=2 ^ ' 

While in ( T4"T|) R could lie either outside or on the surface of the primary, in fj4"8]) it would 
be both convenient and conventional to make R a surface point. Both in (T47I) and (14*8]) . the 

— * 

vector r denotes an exterior point located above the surface point R at a radius r > R 
(with the same latitude and longitude), while r* signifies the position of the tide-raising 
secondary. The quantities fc; are the static Love numbers. 

Under dynamical stressing, the Love numbers turn into operators: 

Utif, t) = (*Y kit) W l {R, f*, t') , (49) 



where integration over the semi-interval t' G ( — 00, t } is implied: 

/D\W pt'=t 

Ui(r,t) = - / Ht-t^WiiR^r*^')^' (50a) 



= f^j [ki{0)Wi{t) - fc,(oo)W,(-oo)] + (^j j' k^t-t') W t (R, r*,t')dt' . (50b) 

Like in the compliance operator (Q]-[2]), here we too obtain the boundary terms: one corre- 
sponding to the instantaneous elastic reaction, ki(0)W(t) , another caused by the pertur- 
bation in the infinite past, — ki(oo)W(— 00) . The latter term can be dropped by setting 
W(— 00) zero, while the former term may be included into the kernel in the same manner 
as in P-E}: 

+ kmwm + (j\ + f k(t-t') w^R^r^t^dt' 

= {j) + f_ ^ - ^(°) + fcj(0)e(*-t')] WtiRr^t^dt' . (50c) 
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All in all, neglecting the unphysical term with Wi{— oo) , and inserting the elastic term into 
the Love number not as ki(0) but as fc;(0) (£ — £') , we arrive at 

Uiftt) = J* k^t-t^W^R, r\t')dt' , (51) 

with ki(t — t') now incorporating the elastic reaction as ki(0) 0(t — t') instead of ki(0) . 
For a perfectly elastic primary, the elastic reaction would be the only term present in the 
expression for ki(t — t') . Then the time derivative of ki would be L(t — t') = k\ 5(t — t'), 
with k\ = ki(0) being the static Love number, and (]5Tj) would reduce to Ui(r, t) = 
( R\ l+1 - 

[yi ki Wi(R, r , t) , like in the static case. 

Similarly to (ITT") , the complex Love numbers are defined as the Fourier images of k l (r) : 

k(x)e lXT dx =k (r) , (52) 

overdot denoting djdr . Following Churkin (1998), the time derivatives k x (t) can be named 
Love functions]^ Inversion of (|52|) renders: 

Hx) = / k(r) z- ixT dr = k t (0) + i X / [*,(r) - fc f (0)e(r)] e'^dr , (53) 
Jo Jo 

where integration from is sufficient, as the future disturbance contributes nothing to the 

present distortion, wherefore fc/(r) vanishes at r < 0. Recall that the time r denotes the 

difference t — t' and thus is reckoned from the present moment t backward into the past 

In the frequency domain, ( )50|) will take the shape of 

/ R \ 1+1 - 

UiU) = [-) Hx) W t (x) , (54) 

X being the frequency, while U l (x) an d W x (x) being the Fourier or Laplace components of 
the potentials U l (t) and W t (t) . The frequency- dependencies k t (%) should be derived from 
the expression for J{x) or p(x) = V^(x) • These expressions follow from the rheological 
model of the medium. 

Rigorously speaking, we ought to assume in expressions (15*21 - "51")) that the spectral com- 
ponents are functions of the tidal mode u and not of the forcing frequency x ■ However, as 
explained in the end of Section HI employment of the positively defined forcing frequencies is 
legitimate, insofar as we do not forget to attach the sign multipliers f[46l) to the terms of the 
Darwin-Kaula expansion for the tidal torque. Therefore here and hereafter we shall expand 
over x i with the said caveat kept in mind. 

5.2 Complex Love numbers as functions of the complex compli- 
ance. The elastic-viscoelastic analogy 

The dependence of the static Love numbers on the static rigidity modulus /x(oo) looks as 

(static) 3 1 

2F1) V7~a\ 



k l = ~1 ~ (sialic) ' ( 55 ) 



12 Churkin (1998) used functions, which he denoted ki(t) and which were, due to a difference in notations, 
the same as our fc ; (r) . 
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where 

(static) (2/ 2 + 4/ + 3) , , 3(2/ 2 + 4/ + 3) , , 3(2/ 2 + 4/ + 3) 



\gpR ^ y ' 4 IttG p 2 R 2 ^ K ' AlnGp 2 R 2 J(oo) 

with p , g, and R being the density, surface gravity, and radius of the body, and G being 
the Newton gravitational constant. The static rigidity modulus and its inverse, the static 
compliance, are denoted here with p(oo) and J(oo) , correspondingly These notations 
imply that we identify static with relaxed. 

Specifically, the static quadrupole Love number will read: 

(static) 3 1 

fc 2 - o ~Tsla~tic) ' ( b '> 

1 1 + A\ 

where the quantity 

(static) _ 57 p(00) 
_ (static) (static) 

is sometimes referred to as p. Clearly, A 2 in (1571) . as well as A h in ( 1551) . is a 
dimensionless measure of strength-dominated versus gravity-dominated behaviour. 

It is not immediately clear whether the same expression interconnects also ki(x) with 
p(x) ■ Fortunately, though, a wonderful theorem called elastic-vis coelastic analogy, also 
known as the correspondence principle, ensures that the viscoelastic operational moduli 
p(x) or J(x) obey the same algebraic relations as the elastic parameters p or J (see, 
e.g., Efroimsky 2012 and references therein). For this reason, the Fourier or Laplace images 
of the viscoelastic equation of motion E and of the constitutive equation look as their static 
counterparts, except that the stress, strain, and potentials get replaced with their Fourier 
or Laplace images, while ki , p , and J get replaced with the Fourier or Laplace images of 

k l (t — t') , P (t — t') , and J (t — t') . For example, the constitutive equation will look like: 
a^ v = 2 p u^ v . Therefore the solution to the problem will retain the mathematical form 
of Ui = kiW[ , with ki keeping the same functional dependence on p , R, and p (or J) as 
in (|56|) . except that now p and J are equipped with overbars: 

Hx) = W=T) i + a)p {x) ,p ^ 



3 1 3 %) 



2(1 - 1) 1 + A l J/J( X ) 2(1 - 1) J( x ) + Ax J 



(59b) 



where 



(2/ 2 + Al + 3)p 3 (2Z 2 + Al + 3)p 3 (2 I 2 + Al + 3) J' 1 
'~ IgpR ~ AIttG P 2 R 2 ~ AlixGp 2 R 2 ' { } 

Although expression f)60p for factors A\_ is very similar to expression (156 p for their static 
counterparts, an important difference between (1561) and (1601) should be pointed out. While 



In the equation of motion, we should neglect the acceleration term and the nonconservative inertial 
forces. Both neglects are justified at realistic frequencies (for details see the Appendix to Efroimsky 2012). 
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in (156]) we had the static (relaxed) rigidity and compliance, /i(oo) and J(oo) = l//i(oo) , 
in (1601) the letters /i and J may stand for any benchmark values satisfying J = 1/fi . This 
freedom stems from the fact that the products Ai J entering (159bl) bear no dependence upon 
J or 11 . The second term in the denominator of (I59al) contains Ji . For convenience, we 
multiply and then divide Ji{x) by some fi , and make the multiplier fi a part of Ai as in 
(1601) . This makes it easier for us to compare (1601) with its static predecessor (156]) . However 
the constant /i in equations (159]) and (160]) is essentially arbitrary, and is not obliged to 
coincide with, say, unrelaxed or relaxed rigidity. Accordingly, J = is not obliged to be 
the unrelaxed or relaxed compliance. 

The above caveat is important, because in certain rheo logical models some of the unre- 
laxed or relaxed moduli may be zero or infinite. This will happen, for example, if we start 
with the Maxwell or Kelvin- Voigt body and perform a transition to a purely viscous medium. 
Fortunately, in realistic rheologies such things do not happen. Hence it will be convenient 
(and possible) to identify the J from (1601) with the unrelaxed compliance J = J(0) emerging 
in the rheological model (1331) . Accordingly, the rigidity fi = 1/J from (I60p will be identified 
with the unrelaxed rigidity /i(0) = 1/J(0) . This convention will play a crucial role down the 
road, when we derive formula (1631 ). 

Writing the / th complex Love number as 



The importance of the products \ki(x)\ sine/(x) lies in the fact that they show up in the 
terms of the Darwin-Kaula expansion of the tidal potential. As a result, it is these products, 
and not k\jQ as some think, which emerge in the expansions for tidal forces and torques, 
and for the dissipation rate. 

In an attempt to preserve the popular notation ki/Q , one may define the inverse quality 
factor as the sine of the lag - see the discussion in subsection 12. 3[ In this case, though, one 
would have to employ the tidal lag ei , and not the lag 5 in the material (which we call 
the "seismic" lag). Accordingly, one will have to write not k\jQ but ki/Qi , where 
1/Qi = sine; . 

Importantly, the functional form of the frequency- dependence sine/(x) is different for 
different I . Thus an attempt of naming sin ei as 1/Q would give birth to a whole array of 
different functions Qi(x) • For a triaxial body, things will become even more complicated - 
see footnote [TT1 To conclude, it is not advisable to denote sine; with 1/Q . 

It should be mentioned that the Darwin-Kaula theory of tides is equally applicable to 
tides in despinning and librating bodies. In all cases, the phase angle q = £i(Xi mpq ) 
parameterises the lag of the appropriate component of the bulge, while the absolute value of 
the complex Love number \k{\ = \ ki(Xi rnpq ) I determines the magnitude of this component. 
The overall bulge being a superposition of these components, its height may vary in time. 

5.3 The tangent of the tidal lag 

In the denominator of (I59aj) the term 1 emerges due to self-gravitation, while A\ J / J(x) = 
Ai\fi{x)\/f I describes how the bulk properties of the medium contribute to deformation 



ki{x) = He[ki(x)] + iZm[ki(x)] = \h{x)\e 
we express the phase lag ei(x) as: 

\ki(x)\ sine/(x) = -Im[ki(x)] ■ 



(61) 



(62) 



18 



and damping. So for a vanishing A\ J/\J(x)\ we end up with the hydrostatic Love numbers 
ki = /; — tt , while the lag becomes nil, as will be seen shortly from (164]) . On the contrary, 



2(1 - 1) 

for very large Ai J/J(x) > we expect to obtain the Love numbers and lags ignorant of the 
shape of the body. 

To see how this works out, combine formulae ( 155]) and (I59b|) . to arrive at 



tane. 



Xm [ki(x)] 



Ai J Xm [J( X 



{Tie [%)] } + {Xm [%)] } 2 + At J Tie [J( X 



(63 



A, 



sm 



a it 



a] 



(z- 1 + z~ a sin ( Of ] 



2 " COS 



( aii\ 



r(i 



2; being the dimensionless frequency defined by (1591) . 

Comparing this expression with expression ( 1551) . over different frequency bands, we shall 
be able to explore how the tidal lag e l relates to the lag in the material 5 (the "seismic 
lag"). 

While expression (163bj) is written for the Andrade model, the preceding formula (I63aj) is 
general and works for an arbitrary linear rheology. 



a 



5.4 The negative imaginary part of the complex Love number 



As we already mentioned above, rheology influences the tidal behaviour of a planet through 
the following sequence of steps. A rheological model postulates the form of J( X ) ■ This 
function, in its turn, determines the form of ki(x) , while the latter defines the frequency 
dependence of the products |fc;(x)| sine(x) which enter the tidal expansions. 

To implement this concatenation, one has to express |fc/(x)| sine(x) via J(x) ■ This can 
be done by combining (1591) with (1621) . It renders: 



\ki(x)\ sine,(x) = -Xm[k t (x)] 



Ai JXm [J( X 



2(1 



{ne[J(x)} +AtJy+ (Xm[J( X )] )' 



a quantity often mis-denoted_j as h/Q . Together, formulae (155)) and (1541) . give us the 
frequency dependencies for the factors |fcz(x)| si ne ;(x) entering the theory of bodily tides. 
For detailed derivation of those dependencies, see Efroimsky (2012). 

As explained in Section HI employment of expressions ( 1621 - [641) needs some care. Since 
both U and k\ are in fact functions not of the forcing frequency x but of the tidal mode 
to , formulae (1621-1641) should be equipped with multipliers sgnu lmpq , when plugged into the 
expression for the Impq component of the tidal torque. With this important caveat in mind, 
and with the subscripts Impq reinstalled, the complete expression will read: 



HXi mpq )\ sine z (x 



At J Xm 



Ax, 



Impq) 



(64) 



Impq J 



2(1-1) 



Tie 



Impq) 



+ At J 



Xm 



J(x 



Impq 1 



2 S g nW i mpg 



14 One can write the left-hand side of (|64[) with ki/Q only if the quality factor is defined through ([26| 
and endowed with the subscript I . 
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a general formula valid for an arbitrary linear rheological model. 

To make use of this and other formulae, it would be instructive to estimate the values 
of A[ for terrestrial objects of different size. In Table [U we present estimates of A 2 for 
Iapetus, Mars, solid Earth, a hypothetical solid superearth having a density and rigidity of 
the solid Earth and a radius equal to 2 terrestrial radii ( R = 2Rq ) , and also a twice larger 
hypothetical superearth ( R = 4_R ffi ) of the same rheology. 



Table 1: Estimates of A 2 for rigid celestial bodies. The values of A 2 are calculated 
using equation ([58]) and are rounded to the second figure. 





radius R 


mean density p 


mean relaxed 
shear rigidity /x(oo) 


the resulting 
estimate for A 2 


Iapetus 


7.4 x 10 5 m 


1.1 x 10 3 kg/m 3 


4.0 x 10 9 Pa 


200 


Mars 


3.4 x 10 6 m 


3.9 x 10 3 kg/m 3 


1.0 x 10 11 Pa 


19 


The Earth 


6.4 x 10 6 m 


5.5 x 10 3 kg/m 3 


0.8 x 10 11 Pa 


2.2 


A hypothetical superearth 
with R = 2 i?^ and the 
same rheology as the Earth 


4.5 x 10 8 m 


5.5 x 10 3 kg/m 3 


0.8 x 10 11 Pa 


0.55 


A hypothetical superearth 
with R = 4 i?^ and the 
same rheology as the Earth 


9.0 x 10 8 m 


5.5 x 10 3 kg/m 3 


0.8 x 10 11 Pa 


0.14 



Taken the uncertainty of structure and the roughness of our estimate, all quantities in 
the table have been rounded to the first decimal. The values of Iapetus' and Mars' rigidity 
were borrowed from Castillo- Rogez et al. (2011) and Johnson et al. (2000), correspondingly. 

In Figure [21 we compare the behaviour of k 2 sine 2 = \k 2 {x)\ sin^x) for the values of A 2 
appropriate to Iapetus, Mars, solid Earth, and hypothetical superearths with R = 2 R^ 
and R = 4 R^ , as given in Table [TJ Self-gravitation pulls the tides down, mitigating 
their magnitude and the value of the tidal torque. Hence, the heavier the body the lower 
the appropriate curve. This rule is observed well at low frequencies (the viscosity-dominated 
range). In the intermediate zone and in the high-frequency band (where inelastic creep 
dominates friction), this rule starts working only for bodies larger than about twice the 
Earth size. If we fix the tidal frequency at a sufficiently high value, we shall see that the 
increase of the size from that of Iapetus to that of Mars and further to that of the Earth 
results in an increase of the intensity of the tidal interaction. For a R = 2 R^ superearth, 
the tidal factor k 2 sin e 2 is about the same as that for the solid Earth, and begins to decrease 
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for larger radii (so the green curve for the larger superearth is located fully below the cyan 
curve for a smaller superearth). 

6 Tidal dissipation versus seismic dissipation, 
in the inelasticity-dominated band 

In this section, we shall address only the higher-frequency band of the spectrum, i.e., the 
range where inelasticity dominates viscosity and the Andrade model is applicable safely. 
Mind though that the Andrade model can embrace also the near-Maxwell behaviour, and 
thus can be applied to the low frequencies, provided we "tune" the dimensionless parameter 
( appropriately - see subsection 13.41 above. 

6.1 Response of a sample of material 

At frequencies higher than some threshold value Xo > dissipation in minerals is mainly due to 

inelasticity rather than to viscosityl^l Hence at these frequencies £ should be of order unity 

or smaller, as can be seen from f!33bj) . This entails two consequences. First, the condition 

X^> 1/ (C t m) > i- e -5 -2 3> 1 is obeyed reliably, for which reason the first term dominates the 

i 

denominator in ( )38l) . Second, either the condition 2 > 1 is stronger than z ^> ( X a or 
the two conditions are about equivalent. Hence the inelastic term dominates the numerator 
in ((38]): z~ a > z~ l C ■ 

Altogether, over the said frequency range, (1381) simplifies to: 

tan5 « ( X r A ) sin (— ) I> + 1) = ( X Cr M ) sin (-) T(a + 1) . (66) 

Clearly tan 5 <C 1 , wherefore tan 5 ~ sin 5 ~ 5 . For the seismic quality factor, we then 
have: 

(seismic)Q-l w (xr j-% in ^ r(a+1) y (67) 

no matter which of the three definitions (12T1- 1231) we accept. Be mindful, that here we use 
the term seismic broadly, applying it also to a sample in a lab. 

6.2 Tidal response of a homogeneous near-spherical body 

Recall that defect unpinning stays effective at frequencies above some threshold Yn , which 
is likely to be above or, at least, not much lower than the inverse Maxwell timeM E.g., 
for the solid Earth, xo ~ 1 J 1 " 1 while r M ~ 500 yr. Over this frequency band, the free 
parameter ( may be of order unity or slightly less than that. (This parameter grows as 

15 For the solid Earth, this threshold is about 1 yr -1 (Karato & Spctzler 1990). Being temperature 
sensitive, the threshold may assume different values for other terrestrial planets. 

Also mind that the transition is not sharp and can extend over a decade or more. 

16 Dislocations may break away from the pinning agents (impurities, nodes, or jogs), or the pinning agents 
themselves may move along with dislocations. These two processes are called "unpinning" , and they go 
easier at low frequencies, as the energy barriers become lower (Karato & Spetzler 1990, section 5.2.3). 
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Figure 2: Negative imaginary part of the complex quadrupole Love number, ki sin €2 = 
— Xm [^(x)] 1 as a function of the tidal frequency x- The black, red, and blue curves refer, 
correspondingly, to Iapetus, Mars, and the solid Earth. The cyan and green curves refer to the two 
hypothetical superearths described in Table [TJ These superearths have the same rheology as the 
solid Earth, but have sizes R = 2i?0 and R = 4R^ ■ Each of these five objects is modeled with 
a homogeneous near-spherical self-gravitating Andrade body with ex = 0.2 and t a = t m = 
In the limit of vanishing tidal frequency Xi t ne factors k-i sin €2 approach zero, which is natural 
from the physical point of view. Indeed, an Impq term in the expansion for tidal torque contains 
the factor ki(Xi mpq ) s ™ e '(Xi mpg )- O n crossing the Impq resonance, where the frequency Xi mpq 
goes through zero, the factor ki(Xi mp q) s ™ e '(^«mpg) mus t vanish, so that the Impq term of the 
torque could change its sign. 
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the frequencies become short of xo •) Under these circumstances, in equation ([S3"]) we have: 
(z -1 <C z~ a <C 1 , whence equation (163]) becomes: 

tane, « YTJi Z ~° ^ (t) 1(1 + a) " (68) 
In combination with (16"B"]) . this renders: 

M 

tane, = — tan 5 . (69) 

1 I + Ai K ' 

Had we defined the quality factors as cotangents, like in ( 123]) . then we would have to conclude 
from (16"91 that the tidal and seismic quality factors coincide for small objects (with Ai ^> 1 ) 
and differ very considerably for superearths (i.e., for Ai -C 1)- Specifically, the so-defined 
quality factor Q t of a superearth would be larger than its seismic counterpart Q by a factor 
of about Aj~ . 

In reality, the quality factors should be used for illustrative purposes only, because practi- 
cal calculations involve the factor \k l (xi mpq )\ s ^ ne i(Xi mpq ) ren dered by ( 165]) . It is this factor 
which enters the Impq term of the Fourier expansion of tides. Insertion of ( 1361 - 1371) into 
( 165]) furnishes the following expression valid in the inelasticity-dominated band: 



\ k i(Xi mpq )\ sine,(x Inw ) 



2(/-l) (A+ 1)5 



X HI being the boundary between the high and intermediate frequencies, i.e., between the 
inelasticity-dominated band and the transitional zone. Expression (170]) resembles the fre- 
quency dependency for |J(x)| sin<5(x) = — 1fn[Jix)] at high frequencies (see equation [361) . 
In Figure 2, dependency ( 1701) corresponds to the slowly descending slope on the far right. 

A detailed derivation of (1701 from formulae (1561 - 157]) and (16"5]) is presented in the Appendix 
to Efroimsky (2012). For terrestrial objects several times smaller than the Earth (so A { ^> 
1 ), the threshold turns out to be 

a 

X H1 = T-jC~ a ■ (71) 
For superearths (i.e., for A x <§C 1 ), the threshold becomes 

X HI = r; 1 = r" 1 C 1 • (72) 

Near the borderline between the inelasticity-dominated band and the transitional zone, the 
parameter ( could be of order unity. It may as well be lower than unity, though not much 
(hardly by an order of magnitude), because too low a value of ( would exclude viscosity 
from the play completely. We however expect viscosity to be noticeable near the transitional 
zone. 

Finally, it should be reiterated that at frequencies lower than some Xo t ne defect- 
unpinning process becomes less effective, so inelasticity becomes less effective than viscosity, 
and the free parameter ( begins to grow. Hence, if the above estimates for x H i turn out 
to be lower than xo , we should set x H i = Xo "by hand" . 
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7 Tidal dissipation versus seismic dissipation, 
in the viscosity-dominated band 

When frequency x becomes short of some Xq i the rate of defect-unpinning-caused inelastic 
dissipation decreases and viscosity begins to take over inelasticity. 

If we simply assume the free parameter ( to be of order unity everywhere, i.e., assume 
that the Maxwell and Andrade timescales are everywhere comparable, then application of 
the Andrade model will set Xo to be of order r^ 1 . Anelasticity will dominate at frequencies 
above that threshold, while below it the role of viscosity will be higher. This approach 
however would be simplistic, because the actual location of the threshold should be derived 
from microphysics and may turn out to differ from r _1 noticeably. For example, in the 
terrestrial mantle the transition takes place at frequencies as high as 1 yr~ l (Karato & 
Spetzler 1990) and may be spread over a decade or more into lower frequencies, as we shall 
see from equation (173|) . 

Another somewhat simplistic option would be to assume that ( ~ 1 at frequencies 
above %o > aric ^ to set ( = oo at the frequencies below Xq ■ The latter would be equivalent 
to claiming that below this threshold the mantle is described by the Maxwell model. In 
reality, here we are just entering a transition zone, where ( increases with the decrease of 
the frequency. While it is clear that in the denominator of (I38p the first term dominates, 
the situation with the numerator is less certain. Only after the condition 

l-a 

C » (Xrj a « ( X r M ) 4 (73) 

is obeyed, the viscous term 1/ (x r M ) becomes leading. This way, although ( begins to grow 
as the frequency decreases below xo > the frequency may need to decrease by another decade 
or more before threshold (1731) is reached. 



whence 



7.1 Response of a sample of material 

Accepting the approximation that the transition zone is narrow^ and that the predominantly 
viscous regime is reached already at Xo or shortly below, we approximate the tangent of the 
lag with 

tan5 « (xr^ 1 , (74) 

(X^/T 1 for r- 1 < X < Xo > 
sin 5 ps { (75) 
1 for 0<x« r- 1 . 

7.2 Tidal response of a homogeneous near-spherical body 

When viscosity dominates inelasticity, expression (163]) gets reduced to the following form: 

tane, ps — ( z~ x = — t — T — , (76) 

1 + A l + (Cz- 1 ) 1 + A i + (xr M ) X Tzif 



17 For a broader transition zone, the rheology will approach that of Maxwell at lower frequencies. This 
though will not influence our main conclusions. 
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comparison whereof with ( I74|) renders: 



M , Ai 

T tano = 

l + A l + { X r M Y 2 l + ^ + tan 2 5 



tane, ~ k- tan § = ^— tan<5 . (77) 



Now two special cases should be considered separately. 
7.2.1 Small bodies and small terrestrial planets 

As illustrated by Table HJ small bodies and small terrestrial planets have A\ ^> 1 . So 
formulae (!76|) and (fT7|) take the form of 



fOT r„ A + 1 ^ * X ° 
Ai x r M for < x < 



tane. 



and 



1 



tane 



r M + 1 t m ^/Ai 



A/ for < x < 



T J" 1UI U \ y ^ , . ~ 1= 

tan<5 r M vC5T+T r M 

Had we defined the quality factors as cotangents of e { and 5 , we would be faced with 
a situation that may at first glance appear embarrassing: in the zero- frequency limit, the 
so-defined tidal Q t would become inversely proportional to the so-defined seismic Q factor. 
This would however correspond well to an obvious physical fact: when the satellite crosses the 
Impq commensurability, the Impq term of the average tidal torque acting on a satellite must 
smoothly pass through nil, together with the Impq tidal mode. (For example, the orbital 
average of the principal tidal torque Impq = 2200 must vanish when the satellite crosses 
the synchronous orbit.) For a more accurate explanation in terms of the |/c;(x)| si ne /(x) 
factors see subsection 17.31 below. 

7.2.2 Superearths 

For superearths, we have Ai 1 , so (177|) becomes 

A[ 11 

tane ; w A iX t m = - — - for 0<x< w — . (80) 

tan 5 r u y /A l + 1 r M 

Here we encounter the same apparent paradox: had we defined the quality factors as cotan- 
gents of e l and 5 , we would end up with a tidal Q z inversely proportional to its seismic 
counterpart Q . A qualitative explanation to this "paradox" is the same as in the subsection 
above, a more accurate elucidation to be given shortly in subsection 17.31 . 

Another seemingly strange feature is that in this case (i.e., for Ai <C 1 ) the tangent of the 
tidal lag skips the range of inverse-frequency behaviour and becomes linear in the frequency 
right below the inverse Maxwell time. This however should not surprise us, because the 
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physically meaningful products ki sine/ still retain a short range over which they demon- 
strate the inverse-frequency behaviour. This can be understood from Figure [2j There, on 
each plot, a short segment to the right of the maximum corresponds to the situation when 
ki sin €i scales as inverse frequency - see formula ( 1ST!) below. 

Thus we once again see that the illustrative capacity of the quality factor is limited. To 
spare ourselves of surprises and "paradoxes" , we should always keep in mind that the actual 
calculations are based on the frequency dependence of \k l (xi mpq )\ s ^ ne i(Xi mp q) ■ 

7.3 Tidal response in terms of |fci(x)| sin e z (x) 

Combining (|3"6"1 - I3"T|) with (|65p . one can demonstrate that in the intermediate- frequency zone 
the tidal factors scale as 

frfr)! sine,( X ) « ^ i)2 {t mX )- 1 , for r~ l » X » (A + l)" 1 , (8b 

which corresponds to the short segment on the right of the maximum on Figure 2. 

From the same formulae ([36] -[37]) and ([65]) . it ensues that the low-frequency behaviour 
looks as 

\k(x)\ sine z ( X ) « 2 {l'-l) A{ Tm X ' fOT T » ( ^ + ^ > X ' (§2) 

a regime illustrated by the slope located on the left of the maximum on Figure 2. 

Details of derivation of (IHTT) and ( [82]) can be found in the Appendix to Efroimsky (2012). 

Just as expression ( [70]) resembled the frequency dependency ( [36]) for | J(x)\ sin5(x) at 
high frequencies, so ([81]) resembles the behaviour of | J(x)\ srn <Kx) a t low frequencies. At 
the same time, ([82]) demonstrates a feature inherent only in tides, and not in the behaviour 
of a sample of material: at x < "^(A + = rj (A + l) -1 ; ^ ne factor si ne Kx) 
becomes linear in X ■ This is not surprising, as the Impq component of the average tidal 
torque or force must pass smoothly through zero and change its sign when the Impq com- 
mensurability is crossed (and the Impq tidal mode goes through zero and changes sign). 



8 Why the Impq component of the tidal torque does 
not scale as R 2l+1 

A Fourier component T lmpq of the tidal torque acting on a perturbed primary is propor- 

21+1 

tional to R k\ sin ei , where R is the primary's mean equatorial radius. Neglect of the 
R— dependence of the tidal factors k\ sine; has long been source of misunderstanding on 
how the torque scales with the radius. 

From formulae ([70]) and ( [HT]) , we see that everywhere except in the closest vicinity of the 
resonance the tidal factors are proportional to A\j{\ + Ai) 2 where A\ ~ R~ 2 according to 
([601 . Thence the overall dependence of the tidal torque upon the radius becomes: 

Over the frequency band x ^> r A ^ 1 (1 + ^) _1 > 

r>2i+i a [ R 2l ~ l , for A t < 1 (superearths) , 
T lm „ ~ an* ~ f— f 2 ~ 

\ ' I R 2l+3 , for A\ ^> 1 (small bodies, small terrestrial planets) 
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In the closest vicinity of the Impq commensurability, i.e., when the tidal frequency Xi mvq 
approaches zero, the tidal factors's behaviour is described by (I82p . This furnishes a different 
scaling law for the torque, and the form of this law is the same for telluric bodies of all sizes: 

Over the frequency band x ^ t - " 1 (1 + Ai)~~ , 

Timpq ~ R 2l+1 h sine; ~ R 2l+1 Ai ~ R 21 ' 1 . I 

9 Conclusions and examples 

Within the inelasticity-dominated band, the phase lags in a homogeneous near-spherical 
body and in a sample of material interrelate as 

j, ( Ai tan S for Ai <C 1 (superearths) , 

tane, = l —- tan 5 ~ < 

1 + A I 

1 [ tan 5 for A[ ^> 1 (small bodies, small terrestrial planets) 

However within the transitional zone, the link between the seismic and tidal dissipation rates 
becomes more complicated. 

The interrelation between the tidal and seismic damping becomes apparently paradoxical 
at low frequencies, where viscosity dominates. As can be seen from (1751 - ISUj) . in the zero- 
frequency limit the tidal and seismic Q s (if defined as cotangents of the appropriate lags) 
become inversely proportional to one another: 



M „ 1 
—— for < X < nr 

tan S r M ^A l 



tane ; « A x xt u = f for 0<x^ 7= • (86) 



This behaviour however has a good qualitative explanation - the average tidal torque Impq 
should vanish on crossing of the Impq resonance. 

While in qualitative discussions it is easier to deal with the quality factors Qi , in practical 
calculations we should rely on the factors k\ sin e/ , which show up in the Darwin-Kaula 
expansion of tides. Just as tan ei , so the quantity ki sin e\ too becomes linear in x f° r 
low values of x- As we saw in subsection I7.3[ this happens over the frequencies below 
X *C 7-t- 1 x . The slight difference between this threshold and the one shown in ( I86I) 

stems from the fact that not only the lag but also the Love number is frequency dependent. 

The factors k\ sin e/ bear dependence upon the radius R of a tidally disturbed primary, 
and the form of this dependence is not always trivial. At low frequencies, this dependence 
follows the intuitively obvious rule that the heavier the body the stronger it mitigates tides 
(and thence the smaller the value of k\ sine; ). However at high frequencies the calculated 
frequency dependence obeys this rule only beginning from sizes about or larger than the 
double size of the Earth, i.e., when self-gravitation clearly plays a larger role in tidal friction 
than the rheology does - see the discussion at the end of subsection I5.4I 

The dependence of k\ sin e\ upon R helps one to write down the overall independence 
of the tidal torque. Contrary to the common belief, the Impq component of the torque does 
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not scale as R 2l+1 , see formulae (1531) and (1541) . 



Here follow some examples illustrating how our machinery applies to various celestial 



1. For small bodies and small terrestrial planets, the effect of self-gravitation is negligible, 
except in the closest vicinity of the zero frequency. Accordingly, for these bodies there 
is no difference between the tidal and seismic dissipations. I 18 l 

Things change in the closest vicinity of the zero frequency. As can be observed from the 
second line of ( 1T5|) . for small bodies and small planets the tangent of the tidal lag be- 
comes linear in the tidal frequency x when the frequency x becomes short of a certain 
threshold!^ x < A « Kr ■ As 

can be seen from ( 1521) . the tidal factor 



ki sine; = srne z(x) becomes linear in x f° r X ^ T M 1 {Ai + 1) 1 ~ r 1 A l 1 . 

2. Tidal dissipation in superearths is much less efficient than in smaller terrestrial planets 
or moons - a circumstance that should reduce considerably the rates of orbit circu- 
larisation. This cautionary point has ramifications also upon the other tidal-dynamic 
timescales (e.g., desp inning, migration). 

In simple words, self-gravity reduces tidal dissipation because gravitational attraction 
pulls the tidal bulge back down, and thus reduces strain in a way similar to material 
strength. 

As can be seen from (180j) . at tidal frequencies x lower than the inverse Maxwell timel^l 
the tangent of the tidal lag changes its behaviour considerably, thereby avoiding di- 
vergence at the zero frequency. According to ( 1521 . the same pertains to the factor 
h sin ei . 

3. While the role of self-gravity is negligible for small planets and is dominant for su- 
perearths, the case of the Earth is intermediate. For our mother planet, the contri- 
bution of self-gravitation into the Love numbers and phase lags is noticeable, though 
probably not leading. Indeed, for [i m 0.8 x 10 11 Pa, one arrives at: 



18 This can be understood also through the following line of reasoning. For small objects, we have A l 3> 1 ; 
so the complex Love numbers (|59bj) may be approximated with 



bodies. 




A 2 



2.2 




3 J(x) 

2 J(X) + A x J 



3 J(X) 



2 Ai J 



+ o(|J/(A ; J)| 2 ) 



The latter entails: 



tane,(x) 



Tie [ki{x)] 



[J(x)] 
Tie [J( X )j 



tan<5(x) 



which is, in fact, correct up to a sign - see the closing paragraph of subsection 15.11 

19 Recall that for small objects Ai ^> 1 . 

20 For superearths, Ai <C 1 . 



28 



so formula (169|) tells us that the Earth's tidal quality factor is a bit larger than its 
seismic counterpart, taken at the same frequency-^ 

(tidal) Q(solid Earth) ~ ^ 5 x (seismic) Q(solid Earth) (88) 

The geodetic measurements of semidiurnal tides, carried out by Ray et al. (2001), 
yield (tidal) q(soiid Earth) w 2 80 . The seismic quality factor ( seismic ) Q( solid Earth ) varies 

over the mantle, assuming values from 100 through 300. Accepting 200 for an arguable 
average, we see that furnishes a satisfactory qualitative estimate. 

This close hit should not of course be accepted too literally, taken the Earth's com- 
plex structure and the uncertainty in our knowledge of the Earth's rigidity. Still, on a 
qualitative level, we may enjoy this proximity with cautious optimism. 



4. The case of the Moon deserves a special attention. Fitting of the LLR data to the 
power scaling law Q ~ \ p has rendered a small negative value of the exponential: 
p = —0.19 (Williams et al. 2001). Further attempts by the JPL team to reprocess 
the data have led to p — —0.07 . According to Williams & Boggs (2009), 

" There is a weak dependence of tidal specific dissipation Q on period. The Q in- 
creases from ~ 30 at a month to ~ 35 at one year. Q for rock is expected to have 
a weak dependence on tidal period, but it is expected to decrease with period rather 
than increase. The frequency dependence of Q deserves further attention and should 
be improved." 

To understand the origin of the small negative value of the power, recall that it emerged 
through fitting of the tidal Q 2 and not of the seismic Q. If the future laser ranging 
confirms these data, this will mean that the principal tide in the Moon is located close 
to the maximum of the inverse tidal quality factor, i.e., close to the maximum taken 
by tane2 in (1751) at the frequency inverse to r M y/Ai . Rigorously speaking, it was of 
course the factor fc 2 sin e 2 which was actually observed. The maximum of this factor 
is attained at the frequency T~ l (A t + l) -1 , as can be seen from (1511 - 1821) . It then 
follows from the LLR data that the corresponding timescale r M {A l + 1) should be of 
order 0.1 year. As explained in Efroimsky (2012), this would set the mean viscosity 
of the Moon as low as 

Vmoou = 3 x 10 16 Pa s , (89) 

which in its turn would imply a very high concentration of the partial melt in the low 
mantle - quite in accordance with the existing models (Nakamura et al. 1974, Weber 
et al. 2011). 

The future LLR programs may be instrumental in resolving this difficult issue. The 
value of the exponential p will have ramifications for the current models of the lunar 
mantle. 

21 When Benjamin et al. (2006) say that, according to their data, the tidal quality factor is slightly lower 
than the seismic one, these authors compare the two Q factors measured at different frequencies. Hence 
their statement is in no contradiction to our conclusions. 
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10 Comparison of our result with that by Goldreich 
(1963) 

A formula coinciding with our (I6"9"]) was obtained, through remarkably economic and elegant 
semi-qualitative reasoning, by Peter Goldreich (1963). 

The starting point in Ibid, was the observation that the peak work performed by the 
second-harmonic disturbing potential should be proportional to this potential taken at the 
primary's surface, multiplied by the maximal surface inequality: 

R R 7 

Epeak ~ ^ 19M i i ~ 19/, + 2g Pj R ' (90) 
2gpR 

R being the primary's radius. 

In the static theory of Love, the surface strain is proportional to R 2 / (19// + 2gpR) . 
The energy loss over a cycle must be proportional to the square of the surface strain. Inte- 
gration over the volume will give an extra multiplier of R 3 , up to a numerical factor: 

R 7 

AE cyde ~ - -2 . (91) 

(19(1 + 2gpRf 

Comparison of (190]) and (191]) rendered 

Q = - ~ (19/, + 2g Pj R) , 

t^E- 1 cycle 

wherefrom Goldreich (1963) deduced that 

Q_ 1 2gpR 

Qo 19 A* ' 

Qo being the value of Q for a body where self-gravitation is negligible. This coincides with 
our formula (16"9"|) . 

In reality, the coincidence of our results is only partial, for two reasons: 

• First, our derivation of the right-hand side of (163]) was based on the prior convention 
that the quantity J entering expression ( 160]) is the unrelaxed compliance J(0) of the 
mantle. Accordingly, the quantity \x = 1/J entering the expression for Ai should be 
the unrelaxed rigidity fi(0) = 1/J(0). In Goldreich (1963) however, the static, i.e., 
relaxed moduli were implied. 

In Ibid., this mismatch was tolerable, because the paper was devoted to small bodies. 
For these objects, Ai is large, no matter whether we plug the relaxed or unrelaxed /x 
into ( I56p . Thence the difference between the tidal and seismic Q factors is small, as 
can be seen from the second line of ( 185]) . 

For earths and superearths, however, the distinction between the unrelaxed and relaxed 
(static) moduli is critical. As can be seen from the first line of (185]) . the tidal Q factor 
is inversely proportional to Ai and, thereby, is inversely proportional to the mantle 
rigidity /i . As well known (e.g., Ricard et al. 2009, Figure 3), the unrelaxed \i of the 
mantle exceeds the relaxed /i by about two orders of magnitude. 
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• Second, as our calculation demonstrates, the simple interrelation given by (1691) and ( )85|) 
works only in the inelasticity- dominated band. In the transition zone (which begins, 
in the solid Earth, at timescales longer than ~ 1 yr) and in the viscosity-dominated 
band of lower frequencies, the interrelation between the tidal and seismic lagging is 
more complicated, and it deviates from Goldreich's formula in a fundamental way. 
In the zero-frequency limit, the cleavage between the tidal and seismic dissipation 
laws gets even larger: the tidal and seismic Q s become not proportional but inversely 
proportional to one another. Description of tidal lagging in all these, low-frequency 
bands requires a rheological model and the subsequent mathematics, and cannot be 
obtained through the simple arguments used by Goldreich (1963). 

Despite these differences, the estimate by Goldreich (1963) provided as close a hit as was 
possible without resorting to heavy mathematics. The elegance of Peter Goldreich's argu- 
ments and the depth of his insight are especially impressive, taken the complexity of the 
problem and the volume of calculations required to obtain the exact answer. 
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Appendix 



Symbol Key 



A t Dimcnsionlcss product emerging in the denominator of the expression for the Love number k t 

E Energy 

E Empirical constant having the dimensions of time, in the generic rheological law l|27|l 

g Surface gravity 

G Newton's gravitational constant 

I Degree (spherical harmonics, Legendre polynomials) 

m Order (spherical harmonics, associated Legendre polynomials) 

J, J(0) Unrelaxed compliance 

J(oo) Relaxed compliance 

J(t — t') Creep-response function (compliance function, kernel of the compliance operator) 

J Compliance operator 

ki Tidal Love number of degree I 

ki(t — t') Kernel of the Love operator of degree I 

ki The Love operator of degree I 

^i(x) Fourier component, at frequency \, of the time derivative of the kernel fc;(t — t ') 

M Mean anomaly 

n Mean motion 

p Exponential in the generic rheological law I I27H 

P t Legendre polynomials of degree I 

Pl m Legendre associated functions (associated Legendre polynomials) of degree I and order m 

Q Dissipation quality Factor 

r Distance 

r Vector connecting the centre of the tidally-perturbed body (interpreted as the primary) with a point exterior to this body 

r* Vector connecting the centre of the tidally-perturbed body (the primary) with a point-like tide-raising secondary 

R Primary's mean radius 

t Time 

u-yv Shear strain tensor 

u-iv (x) Fourier component, at frequency \ , of the shear strain tensor 

U Change in the potential of the tidally-perturbed body (interpreted as the primary) 

W Disturbing potential generated by the tide-raising body (interpreted as the secondary) 

a, /3 Parameters of the Andrade model 

7, v Tensor indices 

r the Gamma function 

8 Material phase lag 
At Time lag 

e Tidal phase lag 

£ lmpq Tidal phase lag of the mode Impq in the Darwin-Kaula expansion 

A Longitude 

£ Parameter of the reformulated Andrade model (ratio of the Andrade timescale t a to the Maxwell time r M ) 

T] Viscosity 

fi, fi(0) Unrelaxed shear modulus (unrelaxed rigidity) 

/i(oo) Relaxed shear modulus (relaxed rigidity) 

fi(t — t ') Stress-relaxation function (kernel of the rigidity operator) 

jl Rigidity operator 

<f> Latitude 

p Mass density 

<r 7 „ Shear stress tensor 

<*7"(x) Fourier component, at frequency \, of the shear stress tensor 

r Time 

t m Maxwell time (viscoelastic timescale) 

t a Andrade time (inelastic timescale) 

T Tidal torque 
0(t — t ') the Heavisidc function 

9 Sidereal angle of the tidally-disturbed body (interpreted as the the primary) 

9 Spin rate of the primary 

f a , <Pu Initial phases of the stress and strain 

X Frequency 

X Frequency threshold marking the boundary between the inelasticity- and viscosity-dominated frequency bands 

Ximpq Physical frequencies of deformation emerging in the tidal theory (absolute values of the tidal modes <*>i m ) 

i^l mpq Tidal modes in the Darwin-Kaula expansion of tides 

u) Argument of the pericentre 

Q Longitude of the node 
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